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Abstract 

A mesoscopic description of spin-transfer effect is proposed, based on the spin-injection mecha- 
nism occurring at the junction with a ferromagnet. The effect of spin- injection is to modify locally, 
in the ferromagnetic configuration space, the density of magnetic moments. The corresponding 
gradient leads to a current-dependent diffusion process of the magnetization. In order to describe 
this effect, the dynamics of the magnetization of a ferromagnetic single domain is reconsidered in 
the framework of the thermokinetic theory of mesoscopic systems. Assuming an Onsager cross- 
coefficient that couples the currents, it is shown that spin-dependent electric transport leads to 
a correction of the Landau-Lifshitz-Gilbert equation of the ferromagnetic order parameter with 
supplementary diffusion terms. The consequence of spin-injection in terms of activation process of 
the ferromagnet is deduced, and the expressions of the effective energy barrier and of the critical 
current are derived. Magnetic fluctuations are calculated: the correction to the fluctuations is 
similar to that predicted for the activation. These predictions are consistent with the measure- 
ments of spin-transfer obtained in the activation regime and for ferromagnetic resonance under 
spin-injection. 
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In the context of spintronics, spin transfer is a generic term that describes the magnetiza- 
tion reversal or magnetization excitations of a ferromagnetic layer provoked by the injection 
of a spin-polarized electric current. This effect has now been observed while measuring the 
magnetization of many different systems [l|]. It has also been implemented into the last 
generation of magnetic random access memories j^j]. 

Spintronics immerged first with the discovery of spin-injection and giant magnetoresis- 
tance (GMR) jjj. Effects of spin-injection at a junction of a ferromagnet are successfully 
described within the two spin-channel model In this description, GMR, or 

spin-accumulation, is due to the spins of the conduction electrons that are driven out-of equi- 
librium by an interface: the difference of the electro-chemical potentials A/i between the two 
spin-channels leads to a redistribution of the spin populations through spin-flip mechanisms. 
Depending on the magnetic state of the ferromagnet at the junction, this redistribution of 
spin-dependent electronic populations modifies the resistance. For this reason, the GMR can 
be used as a probe (e.g. read heads for hard disk drive) in order to measure precisely the 
position of a nanoscopic ferromagnetic moment in the corresponding configuration space. 
Inversely, from the point of view of the ferromagnetic properties, it is natural to expect 
that the spin-flip mechanism at the interface modifies locally the density of the magnetic 
moments in this configuration space. This is a consequence of both the spin-orbit coupling 
(or s-d relaxation), and the redistribution of spins in the magnetic configuration space. The 
spin-injection would then be responsible for a gradient of the density in the magnetic con- 
figuration space, that contributes to the diffusion processes of the ferromagnetic moment. 

The aim of this work is to investigate the consequences of this redistribution mechanism 
for the ferromagnetic moments. This task is performed in the framework of the nonequilib- 
rium theory of mesoscopic systems [9|, |10j. In this context, the coupling between the spin of 
conduction electrons and the ferromagnetic order parameter is due to the introduction of a 
relevant Onsager cross-coefficient that links the spin-polarized electric currents (i.e. trans- 
port of mass, spins, and electric charges in the normal space) to the ferromagnetic current 
(i.e. a massless transport of ferromagnetic moments in the corresponding space). Physically, 
this spin-transfer cross-coefficient accounts for the fact that the spins of the conduction elec- 
trons also contribute to the transport of the ferromagnetic order parameter in the space of 
magnetic moments, in analogy with the thermoelectric Peltier-Seebeck cross-coefficient that 
accounts for the fact that charge carriers contribute also to the transport of heat. 



2 



The motivation for a classical mesoscopic analysis was the need to account for the fol- 
lowing highly specific properties observed in spin-transfer ex per iments that can hardly be 
accounted for in direct microscopic approaches [ll, 
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14j . First, in single domain 



ferromagnets, the reversible part of the hysteresis loop is not significantly modified under 
current injection, while the irreversible jump is drastically modified [141 ] . Second, the am- 
plitude of spin-transfer is proportional to the giant magnetoresistance and to the amplitude 
of the current [3, If], 17]. Third, the Neel-Brown activation law is still valid under current 
injection, with a correction of the barrier height that is quasi-symmetric under both the 



permutation of the magnetic configuration 



change of the direction of the current 
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par allel to anti-parallel and inversely) and the 
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201 ] . "Quasi" -symmetric means here that a 



quantitative shift (a factor 2 to 4 in general) is systematically observed in the amplitude of 
spin transfer for both transitions in spin-valve structures. Fourth, this last quasi-symmetry 
is also observed in the context of ferromagnetic resonance under current injection close to the 
equilibrium states (i.e. with an weak spin- injection) 2l|]. In order to compare the results of 
the model with these observations, after calculating the diffusive correction of the Landau- 
Lifshitz-Gilbert equation (LLG) due to spin-injection in the first section, the correction to 
activation law is derived in a second section, and the correction to the fluctuations is derived 
in the third section. 



A. Spin-injection correction to the ferromagnetic Landau-Lifshitz-Gilbert Equa- 
tion 

In the framework of the two conducting channel approximation, the system is described by 
two electronic populations. A first conducting channel is carrying the conduction electrons 
of spin up t with the electric conductivity (X| and the other channel is carrying the electron 
conduction of spin down j with the electric conductivity <7|. The quantification axis is 
defined by the direction of the magnetization of the ferromagnet. The Ohm's law is valid 
for each channel: in a ID ferromagnetic wire described with a single space coordinate z, the 
electric current (resp. ) is related to the local electric field = — \-j^r Q] through the 
conductivity: J? = cr^i, where /z| (resp. fif) is the electrochemical potential of the channel 
of spins | (resp. I). 

In the case of an interface between a normal metal and a ferromagnet (or between two 



3 



ferromagnets), the spin-dependent relaxation between both electronic populations leads to 
a redistribution of spins within the two channels, that is driven by the difference between 



spin-injection, or spin accumulation, mechanism (see Fig. 1). The equilibrium state of the 
spin system is recovered in the bulk where, by definition, A/z e (oo) = 0. This corresponds to 
a distance of some few times the spin-diffusion length (some tens of nanometers in usual Co 
or Ni layers). The description can be generalized to multichannel model that includes the 
spins of the conduction electron of the s band and of the spins of the conduction electrons 

n □ 

of the d band, with the corresponding interband relaxation [14J, I22j. For convenience, a 
two-channel approximation is used in the following, in which we define Jq = J* + as the 
"spin-neutral" electrical current, and 5J e = J| — Jf is the spin-polarized electrical current. 

The system under consideration is composed not only by the microscopic spins up and 
down carried by the conduction electrons of different nature (band s or d), but also by 
a ferromagnet of length v ( v is also the volume in the normal space for a section unity) 
described by a ferromagnetic order parameter M. The last variable is defined in the space 
of the ferromagnetic moments ( Fig. 2) of constant modulus : M = M s u r , called 7 - space 



in the following [23] . This space can be defined on the unit sphere (see Fig. 2) with the two 
angles 9 and (p, the radial unit vector u r , the azimuth unit vector ug and the zenith unit 
vector u v . The magnetization is then described statistically in the configuration space by 
the density p F (9, (p) of ferromagnetic moments oriented at a given direction 7 = {9, p}, and 
also by the ferromagnetic potential energy V F (9, <p) that contains at least the contributions 
due to the external magnetic field H and the anisotropy energy. Typically, for a uniaxial 
anisotropy with anisotropy constant K, the ferromagnetic potential writes: V F (9, ip) = 
Ksin 2 (9) — M s Hcos(9 — 0) where <fi gives the direction of the applied field. This potential 
energy has the form of a double well potential (Fig. 2). 

In the absence of spin-injection, the magnetization is a conserved variable, so that the 
conservation law writes = —div(j F ), where j F = j F8 ue + j^uip is the ferromagnetic 
current density and the operator div is the divergence defined on the surface of a unit sphere. 
This is no longer the case under spin-injection at an interface: due to the redistribution of 
spins in the different channels (especially from s band to d band) spins are transferred from 



the two electro-chemical potentials A{i e (z) = //? — [l\ (the " spin- neutral" electrochemical 
potential is defined, in turns, by fiQ = //| + /x|). The parameter A/i e accounts for the 
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one sub-system to the other, and the ferromagnetic sub-system becomes an open system. 
In order to work in a larger system that is closed, i.e. that does not exchange magnetic 
moments with the environment, the total ferromagnetic density pf ot and total ferromagnetic 
current j[ ot are defined in what follows. 

In the total system, that includes both the ferromagnetic layer and the spin-polarized 
current, the entropy production d«S/dt (per unit of solid angle and per unit of length) is 
given by 



dS 
~dt 



5,T 



dA/j, e 



J! 



dfJ-o 



(1) 



T \ r ' edz °edz, 
where T is the temperature assumed uniform, V is the gradient defined on the surface of 
the unit sphere, pF is the total ferromagnetic chemical potential, and e the charge of the 
electron. The last term in the right hand side is the Joule heating, the second term is the 
dissipation related to the giant magnetoresistance, and the first term is the ferromagnetic 
dissipation that defines the total ferromagnetic current jf ot in the internal space of magnetic 



moments 



24). 



From the expression of the entropy production Eq. ([I]) and the second law of thermo- 
dynamics dS/dt > 0, the flux involved in the system are related to the generalized forces 
through the matrix of the Onsager transport coefficients 25] 
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All coefficients are known, except the new cross-coefficients / and I, introduced in this 
model, and related to the experimental parameters at the end of the next section. The 



— , and 



electric conductivity o~q is given, in the two channel approximation, by o"o 
the conductivity asymmetry (3 is given by (3 = - ■ The four ferromagnetic transport 
coefficients are {Lee, Le v , L^q, L w }. However, the Onsager- Casimir reciprocity relations 
give Le<p = —L^e, and the symmetry imposes that Lee = L w [14]. Furthermore, the two 
transport coefficients left are not independent since they are both related to the Gilbert 
damping coefficient rj and the gyromagnetic factor V . The following relation holds jl4|: 
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FIG. 1: Illustration of the spin- accumulation occurring in a nanoscopic ferromagnetic layer with 
its two non-ferromagnetic contacts: the profile of the electrochemical potential difference Au e is 
plotted function the spatial coordinate z. 



—aLa 



p F aT 



(3) 



vM s (l + a 2 ) 

where a = 1]TM S is the normalized Gilbert damping coefficient. 

In the absence of spin-injection A// = and there is no coupling between the currents. In 
that case, the well-known LLG equation of the ferromagnetic layer of volume v is recovered 
by inserting the ferromagnetic chemical potential 26] 



p F = k B Tln{p F ) + vV F (4) 

into the expression of the ferromagnetic current density j F = —CVp F where L is the 2x2 
matrix of components {Lee, Lq^ L v e, L^} = Lee/a ({a, —1, 1, a}). We have the expression: 



• F 

JqV 



M s (l + a 2 ) 



U r X VV^ + 



F k B TVp F 



v Po 







j — au r x 


u r x | 



rF + Vp F 

v Po 



(5) 

The LLG equation (that includes diffusion terms) is deduced immediately b y d ividing; 
Eq. (jSJ) with the density p F , thanks to the relations p F ^f = J F and M = M s u r 14]. 
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FIG. 2: The configuration space of the ferromagnetic order parameter is represented by a sphere 
of radius unity. The double well shown is a projection over the plane that contains the two 
equilibrium states at 6\ and 02, and the top of the barrier at 8q. The step approximation for 
the chemical potential fi(9) is plotted with the thermal fluctuations sketched by the dashed area 
around the two minima. 

In the presence of spin injection, A/^ 7^ and a correction to the above LLG equation 
is expected, due to the phenomenological transport spin-transfer coefficient I, introduced in 
the Onsager matrix of Eq. (121) . For the sake of simplicity, we only treat the coupling of 
longitudinal spin relaxation along the vector ug. 

Assuming I constant, the longitudinal component of the total ferromagnetic current in 
the layer J F (6 1 <p, v) = f v j F (9, if, z)dz has the form, 

j eF (e, v , v) = j e /(e, <p)v-ij ^f(e, z)dz (6) 

The quantity 3 Lrg^dz is proportional to the giant magnetoresistance Rqmr generated 



at the interface 



Jj (see Fig. 1 and [29(): 



3 J^jfdz = &eIR GMR {e) (7) 

where Jq = I is the electric current injected in the junction of section unity. Since the 
giant magnetoresistance Rgmr{Q) depends on the angle between the incident spin-polarized 
current (e.g. defined by the magnetization state of a second magnetic layer in a spin-valve 



structure) and the ferromagnetic layer, it depends on the state of the ferromagnetic layer, 
i.e. the position in the 7-space. Eq. fl6]) shows that the total ferromagnetic current is simply 
related to the gradient of a total chemical potential n F ot : 

J F = -£V/if ot (8) 
The total chemical potential jjt F ot writes: 

= + -^(0) (9) 

where the electrospin chemical potential V e (8) is given by an integration over the length 
of the ferromagnet in the direction z, and over the angle 9: 

ye(0) = I f f ^ {Zi e>)dzd6> = [ R GMR {9')de' (10) 
Inserting the ferromagnetic chemical potential Eq. (j3J) into Eq. yields 

= kTln(p F (6)) + -^V e (6)+vV F (6) (11) 

In order to deal with the density of the total system pf ot , the contribution of the spin- 
dependent scattering is included in the logarithm : 

l4* = kTln (pL) + vyF where 

pL(0)=P^0)-e^ Vm (12) 

Consequently, the expression of the ferromagnetic-electrochemical potential takes the same 
form as in the case of a ferromagnetic chemical potential without spin-injection, but with a 
modified density (that depends on the spin dependent scattering process). 
The total ferromagnetic current, given by Eq. (jSJ), writes now: 







| — au r x 


u r x | 



rF , k B TVp 



tot 



V Ptot 



J ~ M s (l + ai)\ UrX [ WV + v p F t 

' (13) 

Once again, the corresponding generalized LLG equation is directly obtained from the 

jF 

expression = = The equation can be rewritten with introducing the "pseudo" effective 

M Ptot 

field H e ff — — VV F + ^f- -jpp* , so that the Generalized LLG equation takes the usual form: 
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du r r r ~ r _ ~ "1 "| 

^ = M a (l + a^) X ^ e// " ^ X K X Heff \ ) (M) 
This equation has the same form than that of the LLG equation without spin-injection 
(note that it is mainly a consequence of the approximation of longitudinal coupling only, 
disregarding precessional coupling), but the diffusion part of the effective field is modified 
through the correction of the density pf ot = p F e k B TL 0s V \ 



, T I Vf>L VaH i dv* 



the correction due to spin-injection o^^^~^ RgmrI is consequently a correction to the 
diffusion term that is proportional to the GMR and to the injected current. However, this 
diffusion term accounts for fluctuations [27]]. The effect of the diffusion terms cannot be taken 
into account in the quasi-static hysteresis loop (i.e. for vanishing temperature or infinite 
measurement times) as a deterministic effective field (this point is discussed in the references 
0, Q])- A deterministic correction to the reversible part of the hysteresis (the quasi-static 
states) is hence not expected here. In contrast, the effect of this correction is considerable in 
the activation regime or in ferromagnetic resonance near equilibrium, described in the next 
sections. 



B. Spin-injection correction to the activation process 

What is the consequence of the diffusion correction Eq. (fT5|) in the activation regime of 
magnetization reversal? This question is investigated below for large time scales (beyond 
nanosecond time scales, or "high barriers"), with the corresponding activation process, i.e. 



the so-called Neel-Brown relaxation [30| . In the activation regime, the effect of the precession 
can be neglected: the gradient V and the divergence operator can be reduced to the scalar 
derivative J|. 

The ferromagnetic potential V F has a double well structure (Fig. 2), with the two minima 
Q\ and 8 2 , and a maximum at 8 . The description of the activation process is based on the 
high barrier approximation under which the ferromagnetic current becomes a step function 
in the 9 space: 
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J F (9, t) = J F {t) [0(9 - d x ) - 0(9 - 9 2 )} (16) 

where 0(9) is the Heaviside step function. 

The chemical potential can also be approximated by a step function that takes the value 
of the equilibrium states in the left (9\) or in the right (#2) side of the potential barrier in the 
ferromagnetic configuration space (Fig. 2) : p F (9,t) = p F (9i,t)O(9 — 9)+fi F (92,t)O(9 — 9 ). 
Using now this expression, the density function p F ot (9) of the ferromagnet in configuration 
space writes: 



pZt(9,t) = P f ot (9i,t)e *S* 0(9, -9)+ fijfi*,t)e ^ 0(9 - 9 ) (17) 

The ferromagnetic current is related to the gradient of the generalized chemical potential 
through our fundamental relation Eq. (JSJ) : J F (9,t) = —Lqq 9 ^"^^ . This equation can be 
written into the more convenient form 

J F (9,t) = -D F ot (9)e~— (18) 
with the diffusion coefficient 

KM = -pT^r = -e hBTLe ° (19) 



where the second equality is deduced from Eq. (jT2l ) and the parameter D F = Lgg-ksT/ p F 
(dimension of angle per unit of time) is the usual ferromagnetic diffusion coefficient that is 
constant. 

The activation process is described by a rate equation (see Eq. (123]) below), e. i. a 
contracted description that is obtained by performing a reduction of the continuous internal 
variable 9 over the equilibrium states d{ (i = {1,2}). 

The total flow has a zero divergence current density divJ F = 0. The system is quasi 
stationary and the total current is X = 2it sin(9) J (t) . Eqs. ([TBI) and Eqs. (ITS]) can be 
integrated over the measure exp(vV/ ksT)d9 to give: 



1 ' 2^m9 [Q{ei ~ 6) ~ Q{ ° 2 -°)\ de = - D L(9) e~^~ ~^~d9 (20) 



Kg J 
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so that the total current writes 



32] 



e fM{e 2 )/k B T _ e n(0i)/k B T 

1 = D ^ 6 ^—j^ 7^ w/kB T ( 21 ) 

h 1 2tt sin(0) aU 

Defining the number of representative points near equilibrium by n(9i) = 
I-^ € p(0')27rsin(6')d6' , the density in the double well potential Eq. (iTTI) leads to the expres- 
sions: 

e M(fli)/fc s r = n to t{0i) 

2ttJ5 e- vVF ( e y^ T stn(e)de 



K02)/k B T _ >Hot\V2) (22 , 



/m^i/At;./' _ n tot (9 2 ) 

" 2tt Jj +e e-v vF ( e )/ k B T sin{6)d6 
where e is a real number. Inserting the above equation into Eq. (12 ID leads to the 
generalized rate equation: 

-7- • //i \ • i a \ n tot{9i) n tot (9 2 ) , . 

X = nto^^i) = -n tot {9 2 ) = — (23) 



Using the steepest descents approximation 



33J for the three integrals present in Eq. ff2T]) 



after inserting Eqs. ( l22l) . the total relaxation times write: 

r M1 = D tot (do)— ^ e * (24) 

More explicitly, this generalized rate equation as the form of the usual Neel -Brown 
relaxation rates T^i±i, with an exponential correction expressed in terms of the electrospin 
chemical potential V e (9i): 

h tot {9 x ) = -n tot (9 2 ) = -^e *b™ m - -^e (25) 

Tl-i-2 T 2 ->1 

where V e (6*) = 8e/ $® RcMR{9')d9' / (3. Accordingly, the Neel-Brown activation law is still 
valid under current injection, with a correction that can be added to the potential energy 
barrier. The total potential energy including the contribution of the spin-injection writes : 

V tot {9) = V F {9) + [ R GMR {9')d9> (26) 

The correction to the energy barrier is proportional to the current I and to the GMR 
integrated over the magnetization states corresponding to the energy barrier height. Eq. 
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(|25|) shows that the process follows the Neel-Brown activation law with the relaxation times: 
r = t e fc s T where AV to t = V to t(0i) — V to t{Qo), Q% = #2} and To is the usual waiting time 
(i.e. the prefactor in Eq. f[2"4"j) ). 

In order to compare this analysis with experimental results performed on nanopillars, 
let us assume a spin-valve structure with two ferromagnetic layers composed of identical 
materials with f3 > 0, in which only two states along the anisotropy axes are allowed. One 
layer is fixed (the pinned layer) and the states of the other (the "free" layer) are investigated. 
The magnetization states of the free layer are Q\ = for the parallel configuration (P) and 
82 = vr for the ant parallel configuration (AP). The GMR for the P state Rgmr{0) = 
corresponds to the reference configuration (no spin-flip). The AP configuration Rgmr{ 71 ') — 
AR corresponds to the maximum GMR. If we take the most simple form for the angular 



dependence of the GMR (3l| R G mr(0) = AR (1 - cos(6)), we have V e (0) - V e (vr/2) 
- 4JA ^- 2) and V e (7r) - V e (vr/2) = + 4/A ^ +2) . 

Note that the GMR parameter AR is a function of f3: expressed in terms of the spin 



diffusion length l s f it writes 14] AR/R = ^rg^f (where v is the length of the layer). 



-P 

In conclusion, the injection of the current leads to suppress one transition and to accel- 
erate the other: the current provokes the magnetization reversal from one configuration to 
the other, and the transition depends on the current direction. In the exemple ebove, the 
current provokes the magnetization reversal from P to AP configuration for positive current, 
and provokes the magnetization reversal from AP to P configuration for negative current. 
This is a sufficient condition in order to accounts for the hysteresis loop of the magnetization 
driven by the current. In the general case, both transitions are defined with the relaxation 
rate 

Ti = T e 2kgT (27) 
with a asymmetry factor c, = RcMR^dO / AR, and the coefficient L e e is defined in Eq. 

( cp = (n — 2)/2 and cap = (vr + 2)/2 in the simple exemple given above). The quasi- 
symmetry under both the permutation of the magnetic configurations and the change of the 
current direction observed experimentally is hence contained in the result expressed in Eq. 

(USD. 

From an empirical point of view, the expression often used in order to fit the data in- 
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troduces the critical current J c (measured at zero external field and extrapolated at zero 
Kelvin), such that: f = tq e - "^ 1- ^ where E a is the anisotropy energy of the ferromagnetic 
layer under consideration. Result Eq. (1251) shows that the critical current is given by the 
expression: 

I,- 



TaE a 



P 



vM s (l + a 2 )d \ el J 8AR 
The Neel-Brown law under current injection is measured in references 
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18, Q. 



(28) 
with 



the typical asymmetry between the two transitions of a factor 2 to 4. The relation I c oc 1/ AR 



is verified in reference 



15[ • The proportionality with /3 is observed through the change of the 



sign while changing the scattering anisotropy 28(. Note that the phenomenological results 
presented here can be generalized to tunnel junctions (see e.g. the work of Schmidt and et 
al. in terms of spin injection in magnetic semiconductors |34|). In the case of tunnel barrier 
a factor ten is typically gained in the magnetoresistance AR, so that the critical currents in 
Eq. fl28l are also decreased by a factor ten [35]. 

What is the value of the spin-transfer coefficient 17 In order to compare with the fer- 
romagnetic transport coefficient Lee/p^ = Ta/(M S (1 + a 2 )) - expressed as the inverse 
of an action J _1 .s _1 - the phenomenological coefficient is compared in the same units: 

l _ E a f3 ( L M \ 
ep^ e8ARI c Ci y p^ J ' 

Experiments performed on typical pseudo spin-valve systems show that it is possible to 
switch the magnetization at zero external field in both directions (AP to P or P to AP in 
the previous exempla) for currents of the order of ± 1 mA. For such currents, the energy 
transferred is of the order of 10 meV 15|: the measured quantity is the slope of the points 
1/ic plotted as a function of AR for the two transitions. The anisotropy energy E a is of the 
order of 0.1 eV. The spin-transfer coefficient l/e is consequently of the order of 10 _1 Lee to 
10- 2 L ee . 

On the other hand, the activation experiments under current injection allow to access 
directly to the spin-transfer coefficient through the Neel-Brown law, without the need to 
measure the activation energy E a . The quantity measured is the slope s = d(ln(r/To))/dI. 



Eq. (EHD and Eq. ([271) show that: ^ = kTsf3 [8e R G¥R 
magnitude 10 _1 is confirmed experimentally in references 
2 to 4 in the asymmetry. 
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-i 



0.1. The order of 



la ], together with the factor 



13 



C. Spin-injection correction to the fluctuations 



Due to its diffusive nature, the correction produced by the spin-injection can hardly be 
observed on the reversible state of the hysteresis. The above section shows that the activated 
irreversible jump of the magnetization is in contrast strongly modified by the spin-injection. 
Beyond the activation process, the presence of supplementary ferromagnetic diffusion pro- 
cesses strongly affects another experimentally accessible parameter: the linear response of 
the ferromagnetic moment to magnetic field, spin-injection, and thermal excitations. The 
response is then proportional to the fluctuations. 

The fluctuations occurring near the quasi-static states in the double-well potential can 



be analyzed from the general fluctuation-dissipation theorem (FDT) in the 7-space 10]. 
The density p F is subjected to random fluctuations that are introduced through a random 
current J* which satisfies FDT: 

(J r F (9,t)J r (9',t')) = 2DUd)(pL(0,t))5(6-6')5(t-t') (29) 
The variation of density is now corrected by the presence of the fluctuation current: 

llPtot{0, t) = -J F (6, t) - J r F (8, t) (30) 



Applying step by step the method described above for the rate equation (Eqs. (120]) to 
to the Eq. (13"U|) , the following expression of the fluctuations is obtained (see reference 



E3D) 

my. 



(X r {t)X r {t')) = {n 9 {9l) ~ {nL){9%) (31) 

where the relaxation times f are that previously defined. This expression has not the 
usual form of a FDT which means that this theorem, strictly valid when fluctuations take 
place around equilibrium states, is not fulfilled. The theorem is restored near an equilibrium 
state (nf ot ) eq = {n F ot )(Qi) or (n F )t ) eg = {n F ot ){62) because transitions from one equilibrium 
state to the other are neglected. We obtain from Eq. ([3TJ) (I r (t)l r (t')) eq = {n ^ q 5(t - if). 
This last expression is valid in the cases of linear ferromagnetic resonance experiments, i.e. 



in a situation where the current is well below the critical current I c defined in Eq. ( T28|) [211 ]. A 
more complicated behavior (highly non-linear) should be expected for strong excitations near 
or beyond the critical current I c in order to interpret the non-linear resonance experiments 

m. 
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In conclusion, in the case of linear ferromagnetic resonance (FMR) measured below J c 
and observed close to one equilibrium state (6\ or 62), a correction to the response of the 
ferromagnet is expected, that takes the same form as that calculated for the activation 
process: 

(IrWrV))" = 2Df (JJy^lZpMl) m _ (32) 

Teq V k Bl ^66 J 

The behavior expected is then surprisingly similar to that predicted in the case of the 
activation process, except that it holds for the amplitude of the linear response and not for 
the transition rates. For (3 > 0, we expect an exponential increase (resp. suppression) of the 
response in the AP state with a positive (resp. negative) current, and an exponential increase 
(resp. suppression) of the response in the P state with a negative (resp. positive) current. 
This highly specific characteristic is in agreement with that observed experimentally in the 
context of FMR measurements under spin-injection below critical current I c , i.e. a situation 
in which the magnetization is close enough to equilibrium states (see results presented in 



reference 



fell)- 



D. Link with microscopic theories 



Before concluding, a last question must be invoked about the relation between the model 
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37|. The 



presented here and the microscopic theories of spin-transfer torque [ll, 
phenomenological transport coefficient I (and of course the known transport coefficient Lgg 
and /3(7o) could formally be defined from the relevant Hamiltonian expression with the help 



of projection-operator formalisms [38[, or any other techniques |39|, |40] that lead to the 
coupled stochastic transport equations of the spin-polarized current and the ferromagnetic 
order parameter in the corresponding configuration space. The difficulty is to manipulate 
on an equal footing a microscopic degree of freedom, the spin of conduction electrons, and 
a collective variable, the magnetic order parameter. This task is far beyond the present 
report, but it is possible to gain some insight with dimension considerations. The physical 
mechanism proposed here for spin-transfer is based on the spin- injection only, that is respon- 
sible for the supplementary diffusion effect of the magnetization through the modification 
of the local densities of magnetic moments (redistribution of spins) . This redistribution of 
spins between the electric sub-system and the magnetic sub-system is governed by specific 
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spin-flip scattering mechanisms (or spin dependent creation-anhilation mechanisms). The 
microscopic approach would define the relevant mechanism and deduce the typical spin- 
transfer relaxation time T tr . The relation between mesoscopic and microscopic approaches 
can consequently be invoked through the relation between the correction of the diffusion 
constant SD[ ot (expressed in dimension of angle per unit of time) and the relaxation time: 
SDf ot oc Sn/rtr, where Sn is the amount of spins transferred from the electric sub-system to 
the ferro magnet. 

In the approach proposed by Berger in a pioneering work, a spin-transfer process 
at the interface is described at the electronic level by a typical spin relaxation time T s d 
calculated from the s-d exchange Hamiltonian. If we assuming that the relevant spin-flip 
relaxation is governed by this mechanism r tr = r sd , we would have 5D[ ot oc l/r sd under 
the relevant hypotheses. In that case, the spin-transfer described here in terms of diffusion 
process would be a consequence, in parallel to GMR effects, of the s-d exchange interaction 
occuring at the microscopic level. 

E. Conclusion 

A description of spin-transfer has been proposed at the mesoscopic level, based uniquely 
on the spin-injection mechanism occuring at the junction with a ferromagnet. The spin- 
accumulation at the interface leads to a local change of the density of magnetic moments in 
the corresponding configuration space. The gradient of density generates a diffusion process 
of the ferromagnetic order parameter, which is responsible for the magnetization reversal. 

The spin-injection has been described at the interface of a ferromagnet by means of the 
usual two-conduction channel model which simplifies the spin-polarized current and the giant 
magnetoresistance analyses. The dynamics of the ferromagnet is treated by means of an out- 
of-equilibrium mesoscopic model with a ferromagnetic current defined in the configuration 
space of uniform magnetic moments. The coupling between the two currents is introduced 
through a new phenomenological Onsager cross-coefficient I that accounts for the fact that 
the spins carried by the electric charges in the normal space contribute also to the transport 
of ferromagnetic moments in the ferromagnetic configuration space. It as been shown that 
the correction to the LLG equation that governs the dynamic of the magnetization comes 
from a diffusion term. We have found that the Neel-Brown activation law is still valid, 
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with a correction to the barrier height that is proportional to the integral of the giant 
magnetoresistances over the ferromagnetic states, from the equilibrium to the top of the 
barrier. The expression of the critical current I c is given as a function of GMR, the damping 
factor and the new spin-transfer cross-coefficient /. Furthermore, the correction to the 
fluctuations is shown to be analogous to that of the activation and is also expressed as an 
exponential term. These results are consistent with the results obtained experimentally in 
quasi-static modes (hysteresis loops as a function of the current), in the activation regime 
under spin-injection, and in linear resonance experiments under spin-injection with / < I c . 
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